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Abstract 



The fluorescence of a single dipole excited by an intense light pulse can 

lead to the generation of another light pulse containing a single photon. The 

influence of the duration and energy of the excitation pulse on the number 

of photons in the fluorescence pulse is studied. The case of a two-level dipole 

with strongly damped coherences is considered. The presence of a metastable 

state leading to shelving is also investigated. 
PACS. 42.50.Dv, 03.67.Dd, 33.50.-j 
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I. INTRODUCTION 

The security of quantum cryptography is based on the fact that each bit of information 
is coded on a single quantum object, namely a single photon. The fundamental impossibil- 
ity of duplicating the complete quantum state of a single particule prevents any potential 
eavesdropper from intercepting the message without the receiver noticing |T[]. In this con- 
text, the realization of a efficient and integrable light source delivering a periodic train of 
pulses containing one and only one photon, would be an important advantage 0. 

The purpose of this paper is to evaluate the reliability of such a source. Assuming that 
a smart eavesdropper can get the information as soon as the number n of photons in the 
pulse is larger than two (see Appendix 1), we define a fractional information leakage fu as : 

k = ^ (1) 

-T n>l 

where P„>i, Pn>2 are respectively the probabilities to get at least one and at least two 
photons. The value of fu has to be close to zero, while the probability Pg = Pn>i to emit 
one photon during the sampling period should be as high as possible; we note that the 
probability to get exactly one photon is (1 — fii)Pe- For a poissonian light source, we have: 



fu = l-{l-P-')Hl-P,) ^ ^ (2) 

It is therefore possible to have a good reliability with an attenuated poissonian light source, 
but Pg has to be very small, which makes the source quite inefficient. A better way to 
have both a good reliability and high emission probability is to design a device with fully 
controlled quantum properties, able to emit truly single photon l^-pT]]. One possibility to 
perform such an emission is to use the fluorescence of a single dipole (e.g. a single molecule 
or a single colored center). As a single dipole cannot emit more than one photon at a time 
leading to antibunching in the photon statistics of the fluorescence light [|T^-[I^ a pulsed 
excitation of the dipole can be expected to produce individual photons on demand 
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In previous works the emitting dipole was generally considered to be a radiatively 
damped two-level system. In the present paper we will rather consider emitting dipoles 
with strongly damped coherences, as it is the case for single molecules or single colored 
center ||11| at room temperature. The decay time of coherences, associated to non radiative 
processes which occur in the picosecond range, is thus much shorter than the population 
decay time, that is typically in the 10 ns range. On the other hand, systems such as molecules 
or coloured centers often have an extra metastable state, which is very long lived and thus 
can induce "shelving" in the emission process. In order to describe these features, we will 
model the emitting dipole using the three level scheme shown in Fig. |l|. 

Owing to the fast damping of coherences, only level populations will be considered, 
and the system's dynamics will be described by using rate equations between the three levels. 
The system can be excited from ground state |1) to excited state |2) with a pumping rate r. 
The decay rate from level |2) to level |1) is F, but the system can also decay to a metastable 
state |3) at rate /5F. The branching ratio + /5) is usually (but not necessarily) very 
small. The emission rate from the metastable state will be neglected {i.e. no photons are 
emitted from level |3)), but we will assume that the system can go back from level |3) to 
level 1 2) with a rate r^. This "deshelving" effect has been observed experimentally and 
may be important under strong pumping conditions. 

The purpose of the present calculation is to evaluate the efficiency of such a system in 
converting a train of classical light pulses into a train of single photon pulses ( "photon gun" ) 
. We will thus assume that this system is excited by a train of light pulses of duration 6T, 
such that T6T <^ 1. The separation between the pulses is denoted by T, with FT ^ 1. For 
ideal efficiency of the source, the dipole should be coupled to a field mode in a microcavity, 
which is then damped to the outside world. Here we will only consider free-space emission 
of the dipole, assuming that the emitted fight is collected by purely passive ways, such as a 
parabolic retroreflector [0]. The corresponding imperfect detection efficiency will be included 
in the present model, but the possible effect of a microcavity will not be considered here. 

In the following section, we will introduce a useful framework for carrying out the cal- 
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culation. Then we will evaluate the quantities of interest, taking into account the detection 
efficiency. Finally we will present numerical results illustrating the behaviour of the system. 



II. THEORETICAL MODEL 



A. General framework 



The evolution of the populations will be described using the diagonal terms of the density 
matrix abb{t,to;a), which denotes the population of level b at time t, starting from level a 
at time to (where a and b may take any value from 1 to 3). For the following it will be 
convenient to define the probability aj^\t,to;a) to go from state \a) at time to to state \b) 
at time t, with the emission of exactly n photons. The quantities a^^, are linked to the 
populations abb by the relation : 

oo 

abb{t, to]a) = Y, (^bb^ (t^ to] a) (3) 

n=0 

The probability density to emit one and only one photon at time t when the system is in 
the state \a) at time to is given by the probability (J^{t, to] a) to be in the excited state at 
time t without any photon emission : 

Pi = r4?(t,to;a) (4) 

The quantities al^\t,to; a) introduced previously can be related through the following re- 
currence relationship: 

ai'^^'\t,to;a) = f ra^\t ,to;a)aS\t,t ■,l)dt' (5) 

Jto 

In other terms, in order to emit (n + 1) photons, the system has to emit the photon (n + 1) 
at time t', and to emit no photon from t' to t. The rate equations for can be written : 

^{t,to;a)^j:r^^a^^\t,to;a) (6) 
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where similar equations hold for dhb, with coefficients Vcb- Using equations 0, ^ and |, it can 
be shown (see Appendix 2) that the rate coefficients r^^^ are related to the coefficients rc& 
by: 

^c6^ = 'T'cb - ^Sbi6c2 (7) 

where 6ai is one if a = i, and zero otherwise. For the three- level system we are considering 
in figure |I|, we thus obtain the following rate equations: 

ail'=r.?^-{l + m4'Ur,ar^ (8) 

4? = -r.ai? + /?r4?-r,4^ (9) 

afi^ = -raff + r,ai? (10) 

The difference between eq. ^ and the original rate equations for populations is the missing 
term proportional to crff in the last equation. This means that the ground level is no more 
filled after the emission of one photon, and ensures the uniqueness of the emitted photon. 
Equations p|-p!0| allow, with the knowledge of the initial state, to derive the different quan- 



tities of interest. We first consider in subsection [I B the ideal situation of perfect collection 



efficiency. Subsection [II C| deals with non-unity collection efficiency, which substantially 



modifies results of subsection [11 B| . Finally we take into account the effect of the metastable 
state in subsection [11 D[ 

B. Two- level approximation with unit quantum efficiency 

We assume first that all the emitted photons are detected (unit quantum efficiency), and 
that the dipole is initially in its ground state Assuming also that /5 <^ 1, we can neglect 
the probability for the system to go to the metastable state in the time interval between two 
excitation light pulses, and set (T33 ~ 0. Eq. §-[T^ therefore reduce to the following system : 

4? = raS? - r4? (11) 
aj? = -r.i? (12) 



whose solutions are, for t < 5T: 

ai'^{t,to;l)^eM-r{t-to)] (13) 

4°^(t,to; 1) = = ^ (exp[-r(t - to)] - exp[-r(t - to)]) (14) 

and for t > 6T: 

aff(t,to;l) = aff(5T,to;l) (15) 

4?(i, to; 1) = exp[-r(i - 5r)]4°^((5r, to; i) (i6) 



The probabihty Pj^^ to emit at least one photon between two pulses (say in interval [0,7]) 
is then: 

pis) ^ r £ a^j {t,0;l)dt ^ 1 - exp{-rST) - exp(-rT)[l - exp((r - r)5T)] (17) 

This probability is of course increasing with the period T, which has to be large compared 
to to assure the emission of the photon (for instance, exp(— FT) = 5 10~^ for a 10 MHz 
pulse train and fa 10ns). We can therefore set: 

pi9) ~ 1 _ exp(-r(5r) (18) 
The probability Pjf^ to emit exactly n photons is given by: 

Pi'^ = E ^i?(r, 0; l)^£dtl^l-r a^j{t\ t; l)dt j rat'\t: 0; 1) (19) 

where the second equality corresponds to the probability to emit the photon n at time t, 
and no photons within [t,T]. In the limit exp(— FT) — > 0, the probability P^^^ is given by 
the following expression, which is well-behaved when r = F : 

P^'^ = (^— ^)2[exp(-F5T) - exp(-r5T)] - exp(-r5r) (20) 
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C. Non-perfect collection efficiency 



In practice, the dipole cannot be separated from the collection system, and the statistics 
of interest is the statistics of the detected events, rather that the one of the emission events. 
Assuming again that the initial state is the ground state, and denoting as t] is the collection 
efficiency [f] = 1 —r]), the probability to collect no photon between [0,T] is : 

oo 

nS^^ = E v'^Pi'^ (21) 

n=0 

Let us introduce the probability daa to reach state \a) without the collection of any photon, 
which is given by : 

oo 

^aa = E ^'^^i^ (22) 
n=0 

From the above definitions, we have no^"* = J2a ^aa- Using a calculation very similar to the 
beginning of this section (see eq. |^ and Appendix 2), the linear differential system for cXaa 
can be shown to be : 

5-22 = = rail - rcr22 + r^crgg (23) 
5"33 = -Fto-ss + PTa22 - rd5-33 (24) 
ail = -rail + ^t^33 + f]Ta22 (25) 

The correction introduced here, compared to equations consists in the addition of a 



term filling the ground state with a rate corresponding to the probability density fjT to 
emit one photon but not to collect it. This term ensures the collection of one and only one 
photon. If the initial state is the ground state, and within the approximations of subsection 
|II B| (/5 <^ 1, (T33 = 0), this system can be rewritten: 

0^22 = rail - rcr22 (26) 
an = -rail + vTa22 (27) 

This system can easily be solved between [0,T], and we find for nQ^"*, in the limit 
exp(-rT) 0: 



nf,^) = an(r/; T, 0; 1) = 77^22(^7; 5T, 0; 1) + an(r/; 5T, 0; 1) 



(28) 



with 



an(r7;5T,0; 1) = ^—^ exp(-r'5r) + ^-^ exp(-r5r) 
a22(r7;5T,0;l) = -^[exp(-r5r) - exp(-r'5r)] 



and 



r' = ^(r + r + 
1 



(29) 
(30) 

(31) 



(r - r)2 + Af]rV) 
r' = ^(r + r - i/(r - r)2 + Af]rT) 

The probabihty IIq^'* allows us to determine the probability 11^^^ = 1 — 11 
at least one photon. It permits also to obtain the probability to collect one and only one 
photon : 



^0^^ to collect 



00 

n=l 



(32) 



We find thus : 



= + ^^^^^n?^) (exp(-r'5T) - exp(-r'5T)) 
rirTST ( r' — w , V — nr , ,„ \ 



,?^(^-p(-r'^T).£;^exp(-/.r)) (33) 

These results correspond of course to the results of subsection |IIB| if = 1. A simpler 
expression can be obtained by considering in first approximation that no more than two 
photons can be emitted during the light excitation pulse. We then have : 



p{9) ^ p{9) ^ pis) ^ 1 _ p{9) ^ p{9) ^ 



(5) _ 1 
2 ~ -L 



(34) 



Equation 21 can then be written as 



J\f = 1 _ Pi9) + r]P^9) ^ ^2^pig) _ p( 



(35) 



and equation ^ as 



,(9) 



ig) 



(36) 
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D. Influence of the metastable state 



In order to study the effect of the metastable state, the 3-level equations given above 
can be solved analytically in the general case, giving lengthy and not very illuminating 
expressions. In physical terms, a short intense pulse will excite the dipole just as previously, 
but now the dipole may end up in the metastable state. Thus the emission of the single 
photon will be delayed by an amount depending on the time spent in the metastable state. 

For definitiveness, we shall consider here the situation where the transition rate jST to 
the metastable level |3) is weak but not completly negligible; this applies in particular to 
single molecules (see next section) . The probability to populate this level when a transition 
occurs form level |2) is -^-^ ~ /?, so the metastable level is reached every i/3Pi)-^ light pulses 
in average. In a way similar to the approximations of subsection ^IBj we can neglect the 



probability to leave the metastable state and to reach it again in the same cycle [0,T]. We 



can therefore neglect the filling term PTa22 in equations 23, and the probability to stay in 



the metastable state in one cycle is exp[— (Fm + rd)T], or for q cycles : 

Pe = exp[-(FM + r,)gT] (37) 

When the system reaches the metastable level, it therefore remains shelved during a mean 
time {Tm + ^d)'^, that will be assumed to be much larger than T. The probability to reach 
this level is approximately /5Pj^^ for each excitation pulse, so the average time it takes for the 
system to find itself in the metastable state is T{pPj^^^)~^ . The number of emitted photons 
is thus decreased by a factor 



T(/3Pi^))(-i) (F + r,)r 

T(/3Pi^Vi + (F + rd)-i /5Pi^) + (r + r,)T 



Even if (3 is small, the factor M can thus induce a reduction of the photon flux. Obviously 
this decrease in the number of emitted photons has different statistical properties than the 
random deletion considered in the previous section [|1^ . One gets now alternatively periods 
where the source is "on", and periods where it is "off". In a practical system, one may 
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consider to use a "deshelving" laser in order to increase r^, and thus to maximize the 
duty cycle of the dipole. 



III. DISCUSSION 

In this section the above model is used to demonstrate the potentiality of a single emitter 
to produce single photons when excited by a light pulse. This potentiality is evaluated by 
the fractional information leakage fn defined in Eq. ([l|). In particular the influence of 
the duration and the energy of the pulse is investigated. The parameters considered in 
the following corresponds to commercially available laser systems for typical emitters such 
as terrylene in p-terphenyl [l^ or Nitrogen- Vacancy colored centers in diamond [|ll| with 
saturation intensity of the order of 1 MW/cm^. Note also that in all the plots discussed 
below the collection efficiency is taken a.s 1] = 0.2, which is a realistic value for an optimized 
passive collection system at room temperature. 

In Fig. 0the ability of the single emitter source to deliver truly single photons is compared 
to an attenuated Poissonian source with the same number of empty pulses. The fractional 
information leakage fu is plotted as a function of the probability Pg of emitting at least one 
photon. The quantity Pg is varied by changing the pulse power while the pulse duration is 
kept constant. When the pulse duration 6T is ten times shorter than the emitter's lifetime, 
it appears that the occurence of pulses with two photons or more is reduced by one order 
of magnitude when a single emitter is used instead of an attenuated Poissonian source. 
Reducing further the pulse duration to 1% of the emitter's lifetime improve the fractional 
information leakage by another factor of 10. 

Fig. shows the influence of the pulse duration ST on the fractional information leakage 
fu for a given excitation peak power {i.e. for a given r). This would correspond to an 
experiment where the pulses are sliced up in a continuous wave laser with fast optical 
modulators. Of course the shorter the pulse the better fu, but when the pulse is too short 
the probability Pg decreases also since the peak power is constant. Note that Pg can exceed 
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the collection efficiency rj for ST large compared to F"^. In this case the fluorescence pulse 
emitted by the dipole contains much more than one photon, so that even after the rj = 0.2 
attenuation the probability of having more than one photon remains larger than rj. 

In Fig. ^ the fractional information leakage fu and the probability Pg of emitting at 
least one photon are plotted versus the pulse power for a given pulse duration. As expected 
short pulses {6T = 0.01/F) require more power to reach a value of Pe around Pe = t] = 0.2, 
since Pg depends only on the pulse energy r6T. But short pulses offer a better fractional 
information leakage fu, owing to the fact that the shorter the pulse, the lower the probability 
of emitting a photon and being reexcited within the same pulse. 

For usual molecules or colored centers, the excited state lifetime is of the order of F~^ = 10 
ns, and typical saturation intensity when focussed on a sub-micron spot are of the order of 
1 mW. For these types of emitting dipoles laser pulses with 6T = 0.1 ns and peak power of 
1 W (ie pulse energy of 0.1 nJ) will already lead to good results. It has to be recalled that 
the incoherent model used here is valid only when the pulse duration remains larger than 
the coherence decay, that is in the picosecond range. This hypothesis prevents the use of 
extremely short and intense pulses, but is fully compatible with the numbers just quoted. 



IV. CONCLUSION 

We have evaluated the efficiency of a single photon source based upon the pulsed exci- 
tation of an individual dipole, in a regime where coherences are strongly damped and thus 
rate equations are relevant. This calculation applies for instance to the excitation of a single 



molecule or a single colored center at room temperature |TT[ . 

With respect to a radiatively damped two-level system, where either an exact vr-pulse 
or a fast adiabatic passage is required 1^-0, the requirement on the pulse intensity is much 
less stringent. This type of system is thus promising for achieving an all solid state single 
photon source operating at room temperature. 

This work is supported by the European IST/FET program "Quantum Information Pro- 
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cessing and Telecommunication", project number 1999-10243 "S4P", and by France Telecom 
- Centre National d'Etude des Telecommunications under the "CTI" project number 99 IB 
784. 



APPENDIX A: LOSS OF INFORMATION OWING TO PULSES CONTAINING 
TWO PHOTONS OR MORE IN A QUANTUM CRYPTOGRAPHIC SCHEME 

We emphasize that fu defined by eq. |l] clearly gives a physically meaningful evaluation 
of the single-photon character of the pulse. We note in particular that when P„>i <^ 1, 
the condition fu < Pi/ 2 is equivalent to the "anticorrelation" criterion a < 1 that was 
introduced in ref. |T^. We give below a few examples that suggest to conjecture heuristically 



that fu also gives a good indication of the information leakage due to the multiphotonic 
character of the light pulses. The quantitative evaluation of fu, which is the main result of 
this paper, obviously does not depend on the arguments given below. 

For the sake of illustration, the information is supposed to be coded in the photon 
polarisation, but the following discussion remains valid for any types of information encoding. 
A simple strategy for Eve to exploit photon pairs is to tap a fraction f/ = 1 — of the beam, 
and to store the corresponding photons. The polarisation of the stored photons is measured 
later on, when Alice and Bob have disclosed the relevant basis information. Assuming that 
the probability to get more than two photons per pulse is negligible, the probability for Eve 
to catch the information is then 27777P2, which is maximum for rj = 0.5 and takes a value 
P2/2. The relative fraction of Bob's information which is known to Eve is then P2/Pn>i, 
which is just fu- For attenuated light pulses, one gets fu ~ P2/P1 ~ Pi/2. The action of 
Eve creates no polarisation errors, and cannot be distinguished from a 50% random loss in 
the transmission between Alice and Bob. 

Another possible, more sophisticated strategy for Eve is to use a fast polarization- 
insensitive quantum non-demolition measurement of the number of photons in each 
pulse, and to defiect every second photon. The polarization of the defiected photons is mea- 
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sured later on, as said above. The fraction of useful bits is thus Pn>i for Bob, P2 for Eve, 
and the information leakage is again fn = P2/Pn>i- This scheme introduces neither polar- 
ization errors, nor apparent loss. It can be nevertheless be detected by Bob if he analyses 
the photon statistics of the light pulses that he receives. 

In presence of high transmission losses between Alice and Bob, for instance in long 
distance or free-space quantum cryptography, both methods can be combined to give even 



more powerful attacks For instance, let us assume that Eve is able to catch the 

light pulses before they go through the transmission line, and to distribute them to Bob 
through her own lossless line. Using the QND set-up. Eve identifies the pulses with more 
than one photon, keep one of them, and redistribute to Bob the remaining photons in order 
to simulate the low efficiency 77^ of the original line between Alice and Bob. In that case, as 
soon as rji < fu, Eve gets essentially all the information and remains undetected. Though 
some countermeasures are possible, it is now clear that attenuated light pulses and high 
transmission losses are a deadly combination for quantum cryptography [ p!9| , p0[1 . 

As a numerical example, when using attenuated light pulses with a typical value Pi = 0.2, 
a fraction fu = 0.1, i.e. at least 10% of the information may leak to Eve. By comparison, 
the single photon source described in this paper will give Pi = 0.1 and fu = 0.002 for ex- 
perimentally reachable operating conditions (TST = 0.01, r = lOOOF, overall efficiency 20%, 
see text for definitions). In the cryptographic situations discussed above, the information 
leakage to Eve is thus reduced by a factor 50 when using the single photon source. 

APPENDIX B: DERIVATION OF THE EQUATION FOR ^33- 

We will derive here the rate equations for the quantity ahb{j]]t,tQ]a), which has been 
defined as : 

00 

ai,b{f}- t,to;a) = J2 V'^^bb^ ^0; a) (Bl) 



n=0 



Using this definition, we have: 
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d_ 

di 



d °° d 



We have, for every n > 0, and from equation H: 



(B2) 



d_ 

di 



^bbiv; t, to; a)|aiJ)(r/; t, to; a) + En=i ^[^2 '\t, to; a)ail\t, t; 1) 
+ llM2~'\t',to;a)§-Al\t,t';l)dt'] 



(B3) 
(B4) 



As we obviously have all\t,t; 1) = 6bi, and using eq. we can rewrite eq. [B3 



d 



-a,,(r/;t,to;a)= t, to; a) + E^Li ^ ['^wF^a \t,to;a) 

+ EcrS^ II r4t'\t\to; a)a^\t,t'; l)dt'] 



(B5) 
(B6) 



Using again eqs. | and eq. ^ becomes 
d 



dt 



abbiv; t, to; a) = 6biiiTa22{t, to; a) + ^ r''^acc{r]; t, to 



(B7) 



which is equivalent to eq. Equation ^ can then be easily obtained by setting r] = 1, since 
abb{l; t, to; a) = abbit, to; a). 
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FIGURES 

FIG. 1. Level scheme. The fluorescence is collected between level | 2) and | 1). Level | 3) is a 
metastable state. 

FIG. 2. Fractional information leakage fu versus the probability Pg of emitting at least one 
photon for rj = 0.2. The dashed and dotted lines correspond to the fluorescence of a single emitter, 
as described in the text, for different excitation pulse durations ST = O.OIF"^ (dashed line) and 
ST = O.OIF^^ (dotted line). The thin solid line is the fractional information leakage fu for a 
Poissonian source. 

FIG. 3. Influence of the pulse duration with r) = 0.2. The pumping rate is kept constant, 
r = lOOF. The thick line is the fractional information leakage fu. The thin line is the probability 
Pe of emitting at least one photon. Both quantities are plotted versus the normalized duration of 
the exciting light pulse VST. 

FIG. 4. Influence of the pulse power. All traces are plotted versus the normalized pumping 
rate r/T for a given pulse duration ST with rj = 0.2. The solid lines correspond to the fractionnal 
information leakage fu, and the dashed lines to the probability Pe of emitting at least one photon. 
These values are given for ST = O.OIF"^ (thick lines) and for ST = O.IF"^ (thin lines). 
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Fig. 2 (Brouri et al) 
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Fig 4 (Brouri et al) 



